On the mean density of complex eigenvalues for an 
ensemble of random matrices with prescribed 

singular values 



Yi Wei 1,2 and Yan V. Fyodorov 2,1 

1 Isaac Newton Institute for Mathematical Sciences, Cambridge, UK 
2 School of Mathematical Sciences, University of Nottingham, UK 



Abstract 



Given any fixed N x N positive semi-definite diagonal matrix G > we derive the explicit 
formula for the density of complex eigenvalues for random matrices A of the form A = UVG 
where the random unitary matrices U are distributed on the group U(N) according to the 
Haar measure. 

1 Introduction 

The question of characterizing the locus of complex eigenvalues for an N x N matrix A 
with prescribed singular values, that is eigenvalues gt > 0, % = 1, . . . , N of A* A, , with A* 
being Hcrmitian conjugate of A, was considered in classical papers by Horn and Weyl [Tj. 
In the present paper we provide a kind of statistical answer to that question. Define G = 
diag(<7i, . . . , <7iv) > 0, multiply the matrix y/G by a general unitary transformation U from 
the left and average over the unitary group U(iV) with the invariant (Haar) measure. This 
construction induces a natural measure on the set of matrices A = U\/G with given singular 
values gi and in this way provides us with the corresponding random matrix ensemble. 
The complex eigenvalues %i of such matrices A will cover generically an annular domain 
Rmin < \z\ < Rmax in the complex plane with some density p(z). In such an approach 
the statistical characterization of the locus of complex eigenvalues amounts to knowing the 
profile of the ensemble averaged value of that density for a given set of singular values. It 
is easy to understand that such mean density can depend only on \z\. The corresponding 
explicit formula is provided in Theorem 2.1 which is the main result of the paper. 

Apart from being a rather nontrivial mathematical problem, understanding the statisti- 
cal properties of complex eigenvalues of the above-mentioned ensemble is motivated by its 
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applications in the domain of Quantum Chaotic Scattering. In this capacity the problem at- 
tracted attention for some time, and a few partial results were obtained previously in several 
limiting cases [21 |3j [4j [6]. Below we give a brief description of the physical context related 
to the problem. 

One of very useful instruments in the analysis of classical Hamiltonian systems with 
chaotic dynamics are the so-called area-preserving chaotic maps, see e.g. [7] and references 
therein. They appear naturally either as a mapping of the Poincare section onto itself, or as 
a result of a " stroboscopic" description of Hamiltonians which are periodic functions of time. 
Their quantum mechanical analogues are unitary operators which act on Hilbert spaces 
of finite large dimension N, and are often referred to as evolution, scattering or Floquct 
operators, depending on the given physical context. Their eigenvalues consist of N points on 
the unit circle (eigenphases). Numerical studies of various classically chaotic systems suggest 
that the eigenphases conform statistically quite accurately the results obtained for unitary 
random matrices of a particular symmetry (Dyson circular ensembles) . 

Let us now imagine that a system represented by a chaotic map ("inner world") is em- 
bedded in a larger physical system (" outer world" ) in such a way that it describes particles 
which can come inside the region of chaotic motion and leave it after some time via M open 
channels. Models of such type appeared in various disguises for example, in (8j OH [10l [11] and 
most recently discussed in much detail in relation to properties of dielectric microresonators 
in [TJ] . A natural mathematical framework allowing to deal efficiently with such a situation 
was suggested in [4], see also [5] and we mention here only its gross features. For a closed 
quantum system characterized by a wavefunction ^ the "stroboscopic" (discrete-time) dy- 
namics amounts to a linear unitary map ^(n) — * ty(n + 1), such that ty(n + 1) = u^(n). 
The unitary evolution operator u describes the "closed" inner state domain decoupled both 
from input and output spaces. Then a coupling that makes the system open must convert 
the evolution operator u to a contractive operator A such that 1 — A* A > 0. It is easy 
to show that one can always choose A — uyl — tt* where the matrix f is a rectangular 
N x M, M < N diagonal with the entries Ty = S^tj ,l<i<N,l<j<M, < tj < 1. 
With u replaced by A, the equation *£(n+ 1) = A^(n) then describes an irreversible decay 
of any initially prepared state \P(0) ^ 0, assuming that external input is absent during the 
subsequent evolution. The complex eigenvalues Zk of the operator A all belong to the interior 
of the unit circle \z\ < 1 and play the role of resonances for the discrete time systems. Let 
us mention that various aspects of resonances associated with quantum chaotic open maps 
recently attracted considerable attention [TJJ EH- 

The relation with the random matrix construction employed in this paper is now obvious. 
It amounts to replacing the true evolution operator u with its random matrix analogue u 
taken from the Circular Unitary Ensemble (CUE) in accordance with the standard ideas of 
Quantum Chaos, and to denoting G = \f\ — tt* . In this way the task of studying resonances 
is reduced to investigating eigenvalues of random matrices A of the specified type. Q 

The first significant random matrix result on such matrices A seems to have appeared in 
[2] where the authors considered so-called "truncations" of random unitary matrices. In our 
notations that case is equivalent to taking g\ = g- 2 = . . . gu = 0, with all the rest N — M 
of gi being equal to unity. In [3] the mean density of complex eigenvalues was derived in 
the limit N — > oo, with M < oo being fixed and all <fc < 1. The work |5J provided some 
general results on the joint probability density of all N complex eigenvalues z%, as well as a 
few formulae for the few-point correlation functions (the so-called marginal distributions) of 
the eigenvalue densities. Those formulae however only lead to explicit managable expressions 



1 Let us however note that the case most frequently encountered in direct physical applications actually 
corresponds to choosing u to be unitary symmetric |12j taken from COE. Such a choice reflects the inherent time- 
reversal invariance typical for the the closed quantum chaotic system. A somewhat simpler choice of unrestricted 
unitary random matrices from CUE corresponds to system with broken time reversal invariance. 
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again in the same limiting case as in [4]. As to the results valid for arbitrary finite N, only 
the simplest particular case gi = g 2 = ■ ■ ■ = Qn-i = 9, 9n < g was so far addressed by a 
variety of methods, see [6] for the most recent account and [3] for an early consideration. 



2 Statement of the Main Results 



Our results show that the ensemble-averaged eigenvalue density function p(z) = ^(\z\ 2 ), i.e. 
indeed depends only on \z\. To write the function ^ explicitly we need to introduce a few 
notations. 

Let s' be the Z-th order elementary symmetric polynomials of g^ i — 1,...,N, e.g. 



!. gl = E<=iSi> s 2 = E 



N 

i<jgigh 



, etc. Let us denote sj. 
Define the following functions of {gi} and the complex variable z: 
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where we used the binomial function 



JV! 



In Eq. 



, we defined a matrix 



Gu] = diag(<7i, . . . , <ft_i, gt+i, ■ ■ ■ ,3jv) and used the dash to denote that the index j can not 
be equal to i in the product. 

The main statement of the paper is that the density of complex eigenvalues can be writ- 
ten in terms of the functions defined above. More precisely we state the following 



Theorem 2.1 Let U <E U(7V) be an element of the unitary group and G — diag((?i, . . . , g^ 
be a fixed positive diagonal matrix, such that 0<gi<---<(?Ar<oo. Let U be distributed 
on U(iV) according to the Haar measure. Then the mean density p(z) of complex eigenvalues 
of the matrix A = UVG is given by 



1 N 

i=l 



(2.4) 



where a = + for \z\ 2 > gi, a = — for \z\ 2 < gi. 



Remark: Exploiting that the function F_ is totally antisymmetric with respect to all gi's, 
we can rewrite the above expression in the form: 



(0 \z\ 2 <gi<g 2 ---<9N 
N 

7f E f a(&) 5i < •■■ <g k < \z\ 2 <g k+ i < ■•• <g N (2.5) 
i=k+l 

gi < .92 • • ■ < 9n < \z\ 2 
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Remark: For N > 3, it is not difficult to show that the eigenvalue density function is 
'smooth' at each \z\ = gt, 1 < i < JV, that is it has a continuous derivative. When N = 3, 
the density function is only 'continuous' at g^ but not 'smooth'. 

In the case of degenerate eigenvalues of the matrix G the density function can be derived 
from Theorem 2.1 by taking the corresponding limits as shown below. 

Corollary 2.2: Suppose the diagonal matrix G has the following degeneracies, 

9ki = ••• = fffci+ii, 9k 2 = "* = 5fc 2 +i2' ■ ■ ■ ' 9k 3 = ■ ■ ■ = 9k s +i s , (2-6) 

which is denoted by the short-hand notation 

G = diag(. . . , [# fel ,. . .,g kl+il ], . . . , [g fes , . . .,g ks+ts ], . . . ) . (2.7) 
Define the following two functions 



JV-l 
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The density function is then given by replacing each Xil=o FAfc+n) in the Theorem by 



F^' and making substitutions Eq. (|2.6|) . 

Proof: Consider the following sum: tp = J2li=o ^A(gk+n) when <7fc = • • • = 9k+i- Taking the 
limit gk+i-i -> 9k+i = 9 yields 

lim ip 

9k + i-1^9k + i—9 

fi— 2 7 \ 



X] F A(5fc+n) 

n=0 

i— 2 



\n=0 
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(<?. 
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g fc+i _i=g fc+i =g 



\g k+ i-i=g k +i=g 



(2.10) 



In the second step we have exploited the linearity of s ; as a function of <?'s. Next, we take 
the limit gu+i-2 — » 9- Repeating the procedure i times, we arrive at 



(2.11) 



Applying the results Eq. (|2.11[) to other degeneracies, we obtain the eigenvalue density func- 
tion for U VG under the condition Eq. (|2.6|) . □ 
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Example 1 (rank-one deviation from unitary matrix): 

In the special case of G = diag(gi, [32, ■ ■ ■ , 9n]), with 32 = • • • = 3n = 3, the above procedure 
leads to an especially simple formula for the mean eigenvalue density: 

(2.12) 



which coincides with the known result [3J [5] . 



Remark: In Eq. (|2.12|) , as 3 — > g\, the density function \P — > 00 on [gi , g] and is zero other- 
wise. On the other hand, the integration of ^> over the region [gi,g] yields one. We conclude 
that in this case the density function is simply S(g— \z\ 2 ), as it must be for a random matrix 
A = U^/g which is simply proportional to CUE matrix. 

Remark: In fact, in our derivation of the main theorem, we can extend the domain of g^s 
to include the origin, i.e. our formula holds for g$ > 0. We illustrate this observation in the 
following example. 

Example 2(truncated unitary matrix): 

Consider the case G = diag([#i, . . .,3a/], [3A/+1, ■••,9n]), where g\ = ■ ■ ■ = g M = and 
3m+i = • • • = 9n = 1- By Corollary, we can write the density function of eigenvalues of 
A = UVG as 

* r A \9i = ---=9M=0,gM+i = ~-=gN = l 

«< 1 -W 2 ) M "'(^) M 1^' <2J3) 

In the first step, we have extended domains of 3 to [0, 00) in Theorem 2.1, and correspond- 
ingly, Corollary 2.2. In fact, when G — diag(0 • Im, In-m), eigenvalues of the matrix A in 
Example 2. coincide with those of (N — M) x (N — M) lower right sub-block of a random 
unitary matrix, also known as the 'truncated' unitary matrix. Same results as Eq. (|2.13p has 
been obtained with a completely different method in [2\. 

Remark: As we can obviously always absorb the 11^(1) phase of N x N complex diagonal 
matrix into U, the domain for the matrix G can be defined on Cj\ The eigenvalue density 
function of UVG, averaged over CUE is then obtained by substituting gi — > \g-i\ into Theo- 
rem 2.1. 



Finally, we compare our formula Eq. (|2.5p with numerical simulations. To this end, we choose 
a fixed diagonal matrix G and generate unitary matrices according to the Haar measure. 
We draw a histogram of the radial part of eigenvalues, \z\, of the matrix Uy/G, see FigQ]To 
compare with the histogram, we define the appropriately modified density function ^ld-^l) = 
2|z|^(|2;| 2 ), which is shown by the solid line. From Fig. [H we observe a very good match 
between our formula Eq. (|2.5p and the results of numerical simulations. 
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Figure 1: Histogram of the radial part of the eigenvalue distribution of 5 x 5 matrices A = Uy/G. 
Here G = diag(l, 4, 9, 16, 25) and U is the Haar-distributed unitary matrix, with sample size 
100000 and bin 0.10. Solid line represents the function Vl/i(|z|) as derived from Eq. (|2.5j) . 



3 Main steps of the proof 
3.1 Colour-flavour transformation 

Our starting expression is the following formula [Ml |T5] for the averaged density of complex 
eigenvalues for a general finite-size non-Hermitian random matrix A 



r ^ 1 v 9 r 9 
p(z) = lim — lim — — ( 

7T k^q dz Zf,— >z dzb 



det 



.det 



k i(z — A) 
i(z-A^) k 

k i(zb — A) 
i(zb — A^) k 



(3.14) 



In our case A — Uy/G, where U € U(iV), G — diag(</i, . . . , gjv) > 0. Averaging over the 
Haar measure on the unitary group U(iV) is denoted by < • • • >u- 

Introduce vectors S a = (s l a ) with complex components and their counterparts \a = (Xa) 
with anti-commuting components (Grassmann variables), for all i = 1, . . . , N and a = 1, 2. 

f s i 

This defines two sets of (graded) vectors ipa = ( ? 

\ Xa 

can be represented as integrals over complex variables s and Grassmann variables x m the 
standard way: 



with a = 1,2. The determinants 



< • • • >t/cx / dU I dS\dS2 exp 



-k(S{Si + S2S2) — i(zbSl S2 + ZbS^Si) 



dxidXidX2d,X2 exp -k(x\xi + x\xi) - i(zx\x2 + zx\xi) 



exp i 



(3.15) 



where we defined dSi dS 2 = ]Y a=l H? =1 ds l a ds' a . Note that by Ea. (l3~T41l . averaging over the 
unitary group in the above expression should be carried out after performing the integral 
over the graded vectors {ip^}. 
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Next we change the order of integration over {tp, tp} and U(iV), which is possible due to 
the fact that U(iV) is compact and the integral is bounded. The integration over the unitary 
group can be performed explicitly by exploiting the colour-flavour transformation discovered 
by Zirnbauer [16) : 



( Qb 


m \ 


Q = 






\ 7 l2 


Qf J 


\ &2 


-If J 



dUexpi^lAijig + iplA}^ = J D(Q,Q) expi + ffWWa \ ■ (3.16) 

Such a transformation trades the integration over U(A), where N can be an arbitrary large 
integer for the integration over a considerably simpler 2x2 graded matrices Q defined as 

(3.17) 

Such Q belongs to a Ricmannian symmetric superspace [T7] of the type AIII| AIII. Here, 77's 
and cr's are anti-commuting Grassmann variables. The so-called boson-boson and fermion- 
fermion blocks of Q are given by 

q b G U(l, 1)/U(1) x U(l) = H 2 and q f G U(2)/U(l) x U(l) = S 2 . (3.18) 

The invariant measure on this domain is defined as 

D(Q, Q) = Sdet JV (l - QQ)dQdQ . (3.19) 

where dQdQ = d 2 qi > d 2 qfdaid<J2dpidp2- 

After the colour-flavour transformation Eq. (|3.16p , we get 



J d^ J 



<■■■>£/ cxF(k) = / ## / £>(Q,Q) exp-(^i,^i) 



where dipdip = dS 2 dS%dxidxidx2dx2 and we defined 




(3.20) 



For a fixed complex number z and a given diagonal matrix G, the integral in Eq. (|3.20|) . defines 
a function F(k) of the variable k. Using the integral representation for Bessel functions we 
can show that F(k) is analytic in the half plane Re k>0. 



3.2 Integration over Q and analytic continuation 

Direct evaluation of the integral over the graded matrix Q followed by the integration over 
%j) in Eq. (|3.20p is very difficult. In fact, it is already a highly involved task in a much simpler 
case where Q is a complex number and ip is a complex vector, see [6] for the corresponding 
calculation in such a case. A natural way out could be changing the order of integration 
in Eq. (l3.20p in order to integrate first over {ip,^} by using the standard Gaussian integral 
formula for graded vectors. However, extra care must be taken in performing such a change. 
To understand this consider the integral involving the boson-boson part of the supermatrix 
Q and the complex vectors Si and S2, 

/bosonic = / dS!dS 2 [ dql e -K(SlS 1 + SlS 2 )-{rz b s\s 2 + lZb SlS 1 )+z qb s\S 1+l q b g t SlSl ^ 
J J\q b \<l 
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where we have omitted the trivial Grassmann integrals. Changing the order of integration 
in Eq. ()3.20|) we arrive at 



?t,|<i 



dS\dS2 e 



-k{S\Si+SIS2) — {iz b S\S2+iz b SlSi)+iq b S\Si+iqbgi S^S'l 



(3.23) 



It is clear that Ibosonic is only well-defined in Re« € (l,oo). For k — > 0, which is the 
limit we have to perform in the very end of the calculation, the integration over the boson- 
boson domain forbids changing integration order in Eq. (|3.20[) . Actually such a problem was 
first noticed in [6], and solved by modifying in a non-trivial way the domain of integration 
Eq. (|3.18p over bosonic variables in the colour-flavour transformation. After such a modifi- 
cation one can actually carry out the required change of integration order for any n > 0. 
We shall however see that one can work in the standard parametrisation Eq. (|3. 18[) in the 
allowed region ReK>l, and then continue to < Rck < 1 exploiting analytic properties of 
the function F(k). 

Let us from now on work in the domain Re k>1. Substituting the transformation 
Eq. (l3.16|) into Eq. ()3.15|) . changing the order of integrations over graded vectors {ip,^} and 
the graded matrix Q and integrating out {ip, ip}, we get 



N 

G(k) = / dQdQ Sdet w (l - QQ) J[ Sdet -1 



z b 
z 



z b 
z 



k - igiQ 



(3.24) 



Performing the integration over the supermatrix Q is still a rather involved technical problem. 
We provide below a few comments related to it. 

The Grassmann variables can be integrated out at any stage, and and we find it con- 
venient to carry out that integration at the very beginning. On the other hand, it turns 
out to be important that the integration over the boson-boson part of Q (i.e. qb) should be 
performed before the fermion-fermion part qf. In fact, a quick inspection of the qf integrals 
in Eq. [|3.24p shows that they diverge logarithmically. However, those logarithmic divergences 
are actually a spurious feature of the colour-flavour transformation. The correct way of treat- 
ing such divergencies when performing any supersymmetric colour-flavour transformation is 
to integrate first over the boson-boson sub-manifold. Then, combining all terms which are 
logarithmically divergent, one can show that the divergent parts cancel each other and the 
result is actually finite. 

To perform the integration over the boson-boson part of Q, we introduce polar coordinates 
qb = Vre w , where re [0,1], 6 £ [0, 2n] so that dq b dq b = drdO. In Eq. (pT2l[|) . G{n) is defined 
for Re k>1. For simplicity, we focus on the real n > 1. And further more, we assume all 
<7j < 1, for i = 1, . . . , N. This is done for convenience only and does not reduce generality 
as we can always scale the G-matrix by the magnitude of the largest eigenvalue, and at the 
end of the calculation to scale it back. Under these assumptions we can integrate over the 
angular variable 9 by residue theorem. In this way the result of the original integration 
naturally splits into a sum of contributions from different residues. 

It is crucial that after integrating over 9, we can show G(k) is analytic in the half plane 
Re k>0. Therefore we are allowed to make the required analytic continuation on G(k) to 
Re k>0. Since both F(k) and G(k) are now analytic functions on Re k>0 and F(k) = G(k) 
on k>1, we conclude that F(k) — G(n). We emphasis that this continuation is possible to 
carry out only after angular integration is performed under the assumption n> 1. 

Knowing that we should take the limit n — > in the end of the calculation, we can use 
the condition k « 1 to simplify significantly the integration over the radial part of qb- It 
turns out that one has to distinguish two essentially different cases: gi > \z\ 2 and gi < \z\ 2 . 
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Each of these two cases yields different result when integrating over which explains why we 
have to distinguish F + {gi) from F_{gi) in the final expression. Integration over the fermion- 
fermion part of Q uses certain properties of elementary symmetric functions but is otherwise 
straightforward. Finally, taking derivatives with respect to Zb and z and letting k — > 0, we 
arrive after straightforward but still cumbersome calculations to the formula Eq. (l2.4p . 



4 Open problems 

In conclusion, we would like to mention a few open problems and possible extensions along 
the lines of the present work. An interesting problem would be to investigate the density of 
complex eigenvalues in the limit N — > oo assuming that the matrix g has a finite limiting 
density v{g) = jj S(g — gi) of eigenvalues gi in an interval of the g— axis. A special variant 
of the problem is to assume that gi are eigenvalues of some random Hermitian matrix with 
rotationally invariant measure. This case is in fact equivalent to the so-called Feinberg-Zee 
problem, see [18], which attracted a considerable interest recently. We hope to address it in 
our future publications. 

Another important extension would be to replace matrices U by unitary symmetric ran- 
dom matrices, or to take them from some other groups (e.g. orthogonal). The corresponding 
colour-flavour transformations are known, but the calculations seem to be extremely chal- 
lenging technically. 
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